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Abstract 

This paper provides a methodology for the central management of risk networks, based on the 
extensive one-dimensional machinery available for one-dimensional risk models. Specifically, we 
define efficient subsidiaries of a central branch, a concept that takes into account the cumulative 
dividends generated by the subsidiaries, as well as the cost of bailing them out by the central 
branch. In the case of deterministic central branches with a single subsidiary, we find closed form 
solutions for the value of the subsidiary to the central branch. Moreover, we extend these results 
in the case of hierarchical networks. In the case of non-deterministic central branches with one 
subsidiary, we compute approximate value functions by applying rational approximations, and 
by using the recently developed matrix scale methodology. 
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1 Introduction 


Motivatiou. In recent years considerable effort has been devoted in the finance and insurance 
literature to model default (or ruin) in the context of a system of interacting firms. The insurance 
giant AIG was almost led to bankruptcy after its subsidiary AIG Financial Products faced con¬ 
tingent claims on its derivatives positions and forced a massive government bailout in 2008. AIG 
Financials Products, liquidated in 2011, will probably remain as a textbook example of a nonviable 
subsidiary whose losses were disproportionate with respect to its prior profitability. 

Similar to the problem of a central branch in need to bail out its subsidiaries is the problem of 
central clearing counterparty, mandated by the Dodd-Frank regulations in response to the recent 
crisis. Just as central branch can be ruined by one subsidiary, a central clearinghouse can default 
if its capital is insufficient to cover the default of one (or several) members. This in itself could be 
a systemic event, so it is a current concern if central counterparties are effectively managed. 

Another example which falls under the umbrella of a central branch is a reinsurance company 
which could face potentially large claims, depending on the primary insurers’ dehcit. 

In this paper we consider a general model of central branch risk networks. A key point is defining 
a notion of efficiency of a subsidiary. Whereas intuitively it is clear that not all subsidiaries represent 
viable businesses from the perspective of the central branch, and some should be closed, it is far 
from clear which criteria should be utilized to achieve that. We will illustrate this by comparing 
several such criteria. 

Let us turn now to the mathematical framework of our study. 

Multi-dimensional risk networks (MRN) are dehned by: 

X{t) = u + ct — S{t) = {Xi{t), t > 0, z e X), 


where I = [1,...,/] is a finite set, the vector u represents the capital of the MRN at time 0, the 
vector c represents a constant cash inflow rate, and S{t) is a process representing cash outflows at 
time t, which may include both Levy and Sparre-Andersen renewal components. 

If no boundary condition is specified, we will call this a free spectrally negative MRN. 

Remark 1. The minus sign comes from the one-dimensional case most studied historically, the 
spectrally negative Cramer-Lundberg process, but the case when X{t) is spectrally positive is also 
interesting. The case of spectrally two-sided X(t) is of course interesting, but harder. 

Remark 2. The simplest case is when the component are i.i.d. compound renewal Sparre-Andersen 
processes, generated by i.i.d. pairs of inter-arrival times and claims = 1,2,... 

Ni{t) 

Xi{t) = Ui + Cit- Si{t), Sift) = ^ 6*]*^ i = 1,....,/, 

1=1 

where Nift) = max{A: : T^'^ := 0 renewal counting processes associated to the 

independent inter-arrivals (see e.g. lAAlO^ ). with intensity Xi = , and the claims sizes 

(i) 

Cj ,j > 1 are nonnegative i.i.d. random variables with arbitrary marginal distribution functions 
denoted by Fi{x),i G I, with finite expectation, denoted m^. However, the restriction to jump 
processes is not essential. 
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Remark 3. After reaching speeial subsets exterior to the state spaee, several continuation/regulation 
mechanisms are possible, like absorption, refleetion, or jumping to the interior. These correspond 
to various possible interactions between the eomponents at times of distress. 

Example 1. A toy example with one absorbing boundary and several reflecting bound¬ 
aries. Consider a central branch which must simultaneously manage several subsidiaries. 

The central branch will keep the subsidiaries solvent by bail-outs until the moment of its bankruptcy, 
or according to a pre-specified rule, e.g. until the subsidiary is deemed non-profitable. The sub¬ 
sidiaries will pay dividends to the central branch. Finally, the expected present (net) value to the 
central branch of a subsidiary consists in the difference between the expected discounted payments 
and bail-out amounts. 

This example suggests the following model: 

Deflnition 1. A central branch (CB) network is formed from: 

1. A unit, ealled central branch, with reserves denoted by Xoft), whose ruin time 

To = inf{t > 0 : Ao(t) < 0} 
causes the ruin of the whole network. 

2. Several subsidiaries Xi{t),i = that must be kept nonnegative or above certain pre¬ 

scribed levels, by transfers from the CB. 

For this network, the boundaries Ui = 0,i = 1,..., I are reflecting and uq = 0 is absorbing. 

Remark 4. There are many applications of the central branch concept: a government/central bank, 
a reinsurance company, an insurance group, a central clearinghouse, etc. 

Another interesting application is that of a coalition or default fund created by several insti¬ 
tutions, for bailing them out when they go bankrupt. Interesting issues here are determining fair 
conditions for merging into (profit participation schemes) and splitting out of the coalition. 

One of the most important activities of central branches is acquiring and closing subsidiaries; 

A critical building block in understanding the management of the CB risk network is modeling 
one subsidiary managed by the CB. The benefit is represented by the dividends received from the 
subsidiary, while the costs are the cumulated bail-out costs, the number of bailouts, up to a certain 
time set by the CB and referred to as “patience”. Establishing decoupled rules for the management 
of subsidiaries is the key step towards solving problems pertaining to the whole risk network, such 
as maximizing the total cumulated net income from all the the subsidiaries, until the eventual 
bankruptcy of the CB. 

Notation. Denote the ruin times and rnin probabilities (finite time and eventual) of the 
components when isolated from the network by 

Ti{ui) = inf{t > 0 : Xfit) < 0}, i = 1, 2,... 

= P{Ti{ui) < t), := P{Ti{ui) < oo). 

The ruin probability of the CB and its Laplace transform will be denoted respectively by 

^{t,u) = T(Lri,c) = Pu[to < t], fiq{u) := fiq{u,c) = Eu . 
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One-dimensional first passage problems have been very extensively studied; typically, 
Laplace transforms are available, especially when either Ai or Cj have a matrix exponential distri¬ 
bution, and explicit inversion of the Laplace transforms is also possible sometimes, especially when 
at least one of Ai or Ci have an exponential distribution. 

The eight pillars of one-dimensional first passage problems for spectrally negative 
(and positive) processes. In the last decade it was realized that the solutions of a large gamut 
of problems (dividends, drawdowns, barrier options, exotic options, etc...) may be expressed in 
terms of the solutions of eight basic problems (130]) . ([M]), p4l) . (fill) . (1381) . (j42]), flHO]) . (1^ . which 
at their turn may be ergonomically expressed in terms of a couple of basic solutions called scale 
functions |§ . 

These eight pillars expand the content of the “basic dictionary entry” for a specific Markovian 
stochastic process furnished by its generator, or by its equivalents: the symbol and the first scale 
function, and act as a sort of a instruction kit for solving a variety of first passage problems. 

The precursor of this idea may be found in [AKP04j . where the ’’second scale function” was 
introduced. Its full development may be found in [PisOhj . where a list of identities expressible 
in terms of the two scale functions was offered -see also |KKR13j . Remarkable extensions which 
involve “relaxed” first passage times were provided in [AI14t IAIZ14j . It seemed probable from the 
start that these identities could be applied to Markov processes with state dependent generators as 
well, but the difficulty to compute the scale functions stopped further research, with the notable 
exception of jACP'*~l()) . 

The idea of constructing solutions out of basic scale functions came finally to full bloom while 
being extended to spectrally negative Markov additive processes (which include Levy processes in 
a modulated Markovian environment, and the case of Sparre-Andersen processes with phase-type 
arrivals) in |KP08t Uvallt IIP12] . reviewed in the appendix. This extension requires matrix scale 
functions, and the end result is a Mathematica matrix scale package offered by |Iva m- 

From one to several dimensions. Multi-dimensional first passage problems are considerably 
harder than one dimensional ones, and one cannot expect general formulas [§ . In this paper we 
provide approximations in the multidimensional case by reduction to one-dimensional results. 

Efficient subsidiaries. A crucial issue for a coalition is how to accept efficient members and 
eventually reject them if they are not efficient. For that it is natural to evaluate each member 
separately, by classic one dimensional risk measures like ruin probabilities, or the value of future 
dividend payments made to the coalition. However, the choice of an economic principle for evalu¬ 
ating efficiency is not at all obvious. 

One often used approach is to optimize discounted dividends of rate d = c'j, 7 € (0,1] taken 
above a constant threshold b (7 represents the proportion of income taken above the threshold) 
- see for example [AHT071 IAva09] . The process obtained by subtracting the dividends is called 
reflected when 7 = 1 , and refracted when 7 £ ( 0 , 1 ). 

Judging efficiency as readiness to pay dividends. Another method to judge efficiency 
was made in the conference paper |AM15j . who proposes to define efficiency as local optimality 
of 6 = 0 over some interval [0, e), e > 0. 

The motivation is that such subsidiaries are functional from the start and can contribute cash- 

Mhe eight pillars are respectively the solutions of the one, two-sided and drawdown constrained smooth first 
passage problems, the formulas of reflected bailouts, dividends, severity of ruin and regulated severity of ruin, and 
the resolvent density. 

^One exception is a Pollaczek-Khinchine type formula for the transform of ruin probabilities V['(ii) of spectrally 
negative networks provided in the foundational paper | CYZ03| . However, this formula involves several unknown 
functions (the Laplace transforms over each boundary facet of the state space), and it isn’t at all obvious how to 
exploit this formula numerically. 
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flows to the central branch without having to wait first until its reserves build out; effectiveness 
is thus translated in this paper as readiness. A second motivation is that inefficient subsidiaries 
may be turned into efficient by setting up rules to monitor their time passed in in ” orange zones”, 
and eventually close them when necessary. This may be achieved by “killing” them with a rate 6i 
in the orange zone, where 9i is chosen to render the barrier 6 = 0 locally optimal. Furthermore, 
one may use ’’two step” killing rates, and choose the killing rate as oo in ’’red zones”, as suggested 
in [CR15] . We show in a variety of example that this procedure produces reasonable results. 

Summing up, we will propose that a subsidiary is: 

1. Non-efRcient and rejected immediately if its loading factor ~ 1 •= P~^ ~ ^ is not 

nonnegative, since this implies an infinite number of bail-outs. 

2. Totally efficient and accepted for ever in the coalition iff 

p = \E[Ci]/c < 1, and k < f{q) (1) 

where q is the discount rate, and f{q) is an increasing function of q, with f{q) < 0, obtained 
as optimality of 6 = 0 for some specihc dividends distribution scheme, and A: > 1 captures 
the cost associated with capital infusions towards a subsidiary. 

3. Partially efficient if the loading condition p < 1 is satisfied, and k > f{q). These subsidiaries 
will also be accepted, but only with an impatience rate 9, resulting in killing the subsidiary 
after its time or lowest value in an orange zone exceeds an exponential r.v. of rate 9. The 
impatience rate 9 is chosen so that 

f{6,9) = k, 

where f{5,9) is computed from the optimality of 6 = 0 for the impatience modified value 
function. This is illustrated in example [3] below, where f{6, 9) = f{q + 9). 

The classic De Finetti objective maximizes expected discounted dividends until the ruin 
time. For judging efficiency, it is natural to take also into account a final bail-out, resulting in the 
optimization objective: 


V^^\x) 


SUpFj; 

TT 



e-'i^dD^it)+ e-'i^w{U{T)) , 


( 2 ) 


where w{u) is a so called Gerber-Shiu penalty function. 

The optimal dividend distribution is of multi-barrier type [Ger72] . and the end result may be 
expressed in terms of scale functions [APP071 IAPP15| . Further conditions are necessary to ensure 
that single constant barrier strategies suffice |AM05[ [APP071 fLoeOSl ILRIO] . 

Judging efficiency of subsidiaries over an infinite horizon by optimizing bail-outs 
and dividends. Over an infinite horizon, the subsidiary will possibly need to be bailed out a 
number Nb > 0 of times. 

In the case of linear transaction costs ku — K, the optimization objective (of particular interest 
in a bail-out setting) becomes the expectation over an infinite horizon of a linear combination of 
discounted dividends D{t), cumulative bailouts Z{t), and number of bailouts N'^{t) up to time t: 


V^^\x) 


sup Fa; 

TT 



e-^^dNl{t) 


(3) 


Since in a diffusion setting this objective has first been considered by Shreve, Lehoczky, and 
Gaver (SLG) |SLG84] - see also Lokka and Zervos [TZOSj - we will call it the SLG objective. 
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For spectrally negative Levy processes, the optimal dividend distribution for the SLG objective 
is always of constant barrier type, and the end result may be expressed in terms of scale functions 
[APPOTj . 

Bail-out intervention times. For bail-out times, one may consider the classic ruin time 
r = T(j”, and also several interesting alternatives generalizing it: 

1. One may replace r in Q by a ruin time observed with Poissonian frequency r- see for 
example |AIZ14[ IAI15] , which is equivalent to a Parisian ruin time with exponential grace 
period - see for example |LRZ13| . 

2. One may further replace r in ([9]) by a bankruptcy time Tr^a involving Poissonian-Parisian grace 
period below 0 and also absolute killing at a < 0 - see for example |Renl4[ Sec. 4]. 

Some definitions of efficiency based on prior literature. 


Example 2. A simple, but unsatisfactory definition of efficiency. Instead of our notion of 
readiness, consider defining efficiency as nonnegativity of the SLG objective for the 0 barrier, when 
dividends are D{t) = dt. By fAPPOll Thm. 1, ( 4 - 4 )] 

y(0(o) = >0 (4) 

yields k < = p~^. This generalizes easily for the 0 threshold, yielding 


k < 


d 

Ami 



7- 


Example 3. SLG readiness. Consider now efficiency defined as the optimality of b = 0 for 
the SLG objective with constant reflecting barrier. This problem is fully analyzed in \APP01 , 
Thm. 3], and in particular \APP01/\ Lem. 2] shows that the optimal SLG constant barrier is b* = 0 
^ffk<l + ^. 

By this criterion, a subsidiary i with 


ki<l + ^ (5) 

will be deemed totally efficient and accepted for ever in the coalition [§ . Subsidiaries with loading 
condition Pi < I and 

ki > 1 — 

will be deemed partially efficient and accepted only for a random time with law £{0i), where 


9i -\- q — Xi{ki — 1) 


(7) 


(rendering thus b* = 0 optimal with respect to the total discount rate 6i = q + 9i). 

Unfortunately, the criterion © does not take into account the claim size law. 

^Otherwise |APP071 (5.6)] 

b* = inf{a > 0 : [kZ^'^fa) - - kqW‘-'^\af < 0} > 0 (6) 

(i.e. the subsidiary must be forgiven from paying dividends in between (0,6*). See also |KS08I [ES11 | for related 
results. 
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To avoid these shortcomings, we introduce a new efficiency concept in Section 13.11 
Contents and contributions. Our first contribution is to provide in Section [3] a definition of 
efficiency which is acceptable for a wide variety of subsidiaries satisfying the conditions of Lemma 
1, which include exponential claims. 

The central branch model is described in more detail in Section [H Our second contribution, 
Theorem [2] in Section [U applies to a purely deterministic CB Xo(t) = uq + cot with one 
subsidiary. In this case, the computation of the finite time ruin probabilities and other performance 
measures (including total subsidiary dividends until ruin) reduces to the corresponding computation 
for a subsidiary with modified initial capital uo + ui/ki and initial income rate cq + ci /fci . More 
precisely. 


T(t,it,c) = Ti(t,no/A; + ui,co/k + ci), 

V{u,c) = Vi{uo/k + ui,co/k + ci), (8) 

without any distributional assumptions! For example, the subsidiaries may be dual risk 
processes, or spectrally two-sided Levy processes,... 

The proof of this result, via a pathwise argument, yields also an upper bound when / > 1, and 
an extension to hierarchical networks is given in Corollary [TJ The notion of hierarchical networks 
has been used ever since Gerber ([?]) to model chains of reinsurers and is the most important 
application of the deterministic CB case. 

Our third contribution, in Section [6Tl deals with non-deterministic CB’s with one subsidiary, 
which do not admit an exact solution due to their complex dependent Sparre-Andersen structure. 
We isolate the CB and one subsidiary, and allow for a non-deterministic component of the CB 
process to represent other claims of the CB, for example aggregate net flows of the other subsidiaries 
or other liquidity needs of the CB. We propose an approximation approach based on the classic 
idea of rational approximations to replace the Sparre-Andersen structure by a Markov modulated 
Levy structure, and by using subsequently the matrix scale methodology. 

We propose to obtain SNMAP (spectrally negative Markov additive process) approximations 
for non-MAP central branches, by using bivariate phase-type approximations for the joint law 
of the downward ladder time and height. The advantage of this approach is that once a 
SNMAP approximation is obtained, many similar problems may be solved just by applying the 
scale matrix methodology developed by |KP08l llvalll IIP12] . and using the SNMAP Mathematica 
package of J. Ivanovs [Ival3| . Different problems are thus solved simultaneously! 

A numeric illustration is performed in section 16.21 where we consider the problem of choosing 
a barrier B maximizing CB dividends until ruin, in the case of one subsidiary with exponential 
claims (and without dividends). 

In this case, univariate phase-type approximations of the downward ladder density, obtained 
via a continued fraction expansion -see section lO.ll H , provide a SNMAP approximation of the CB 
process. Subsequently, using the SNMAP package of Ivanovs provides the optimal barrier. 

Our methodology based on the three ideas proposed above lays the grounds for an approxi¬ 
mate optimization of risk networks. The first step consists of setting impatience parameters 6i 
(but not the initial allocation capital infusion from the CB Vi and dividend barrier bi) for each 
subsidiary viewed in isolation from the network. The second step will consist in setting dividend 
barrier parameters bi and optimal allocation parameters Vi, using a decoupled objective of the form 
V{v) = Vi{ui + Vi), representing the sum of benefits to the network from all components, under 

^For other potentially useful explicit computations see |APP08I IBCRll) . 

^Note that even though here the CB claim arrival times have an explicit Bessel-type law, our approach replaces 
this by a phase-type approximation. 


7 













the constraint = '^Oj keeping the 6i, and using the barriers 6* = 0 from the hrst step as 

initial values for an iterative procedure. When transaction costs are present, they may also be 
incorporated via the reduction result ([8|). A harder problem would be to take into account the 
possible bankruptcy of the central branch. 

Other alternatives of efficiency based on threshold strategies, linear reflecting barriers, Poisson 
evaluation periods, and killing based on total current bail-out will be investigated in future work. 
To prepare this, a review of the SNMAP matrix scale approach is provided in the last three sections. 


2 Preliminaries on first passage 

Consider a spectrally-negative MAP {X{t),J{t)) and assume that none of the underlying Levy 
processes XI is a.s. non-increasing (see Section [7] for a review of these processes). Dehne the hrst 
passage times by 

= inf{t >0:X{t)>b}, '^b — hT{t > 0 : A(t) < 6}, 

with inf 0 = -|-oo. We will sometimes write r for the ’’ruin time” rj". 

Let Ta b = Ta^b = ^'^b~ denote the “two-sided” exit time from an interval [a, b]. 


If the process is only observed at the arrival times Tr = {Ti,i = 1,2,..., of an independent 
Poisson process of rate r, the analog concepts are the stopping times 


T+ = inf{r, : X{T,) > b}, = inf{t > 0 : A(T,) < 6} 

see for example [AIZ141 IAI15| . 

A related concept is the Parisian ruin time Tr{a) dehned as the hrst time when the most 
recent excursion below a has exceeded an exponential rv Sr of rate r (in other words, the process 
is killed below a at rate r) - see for example [LR.Z13] . When a = 0, the notation will be used. 

The classic De Finetti objective maximizing expected discounted dividends until the ruin 
time may be taken as a basic principle for judging efficiency. In this context, it is natural to take 
also into account a hnal bail-out, resulting in the optimization objective: 




supFj, 

TT 







(9) 


where w{u) is a so called Gerber-Shiu penalty function. 

The optimal dividend distribution is of multi-barrier type [Ger72] . and the end result may be 
expressed in terms of scale functions [APPO?) IAPP15| . Further conditions are necessary to ensure 
that single constant barrier strategies suffice |AM05[ IAPP071 fLoeOSl ILRIO] . 

Over an inhnite horizon, the subsidiary will possibly need to be bailed out a number Nb > 0 
of times. 

In the case of linear transaction costs ku — K, the optimization objective (of particular interest 
in a bail-out setting) becomes the expectation over an inhnite horizon of a linear combination of 
discounted dividends D{t), cumulative bailouts Z{t), and number of bailouts N'^{t) up to time t: 


V^^\x) 


sup Fa; 

TT 



e-^^dNl{t) 


( 10 ) 


Since in a diffusion setting this objective has hrst been considered by Shreve, Lehoczky, and 
Gaver (SLG) |SLG84j - see also Lokka and Zervos [LZn8| - we will call it the SLG objective. 

























For spectrally negative Levy processes, the optimal dividend distribution for the SLG objective 
is always of constant barrier type, and the end result may be expressed in terms of scale functions 
[APPOTj . 

Bail-out intervention times. For bail-out times, one may consider the classic ruin time 
r = T(j”, and also its ’’soft” Poissonian-Parisian alternatives generalizations. 


3 Efficient subsidiaries 


As argued in the introduction, we will consider that a company is: 

1. Totally efficient and accepted for ever in the coalition if 

p = \E[Ci]/c < 1, and k < f{q) (11) 

where q is the discount rate, and f{q) is an increasing function of q, with f{q) < 0, obtained 
as optimality of 6 = 0 for some specific dividends distribution scheme, and A; > 1 captures 
the cost associated with capital infusions towards a subsidiary. 

2. Partially efficient if the loading condition p < 1 is satished, and k > f{q). These subsidiaries 
will also be accepted, but only with an impatience rate r, resulting in killing the subsidiary 
after its time or lowest value in an orange zone exceeds an exponential random variable of 
rate r. The impatience rate r is chosen so that 

f{q,r) = k, 

where f{q, r) is computed from the optimality of 5 = 0 for the impatience modified value 
function. In example [3] we have f{q) = 1 ^ (see ([5])) and we set f{q,r) = f{q + r). 

Motivated by the fact that the criterion ([5]) does not take into account the claim size law, we 
introduce in Section [Q an efficiency concept which does not have this shortcoming. 

We list now several alternatives of dividend payment strategies of a subsidiary towards a coali¬ 
tion, which could turn out to provide appropriate definitions of readiness in our context, and which 
we plan to investigate in the future. 


1. Linear reflecting barrier strategies consist in specifying a function b{t) = b+ {c — d)t, and 
taking as dividends all surpluses above it. 

2. Two-step premium/refraction/threshold dividend strategies postulate that only a pro¬ 
portion 0<7 = ^<lof the premium income may be paid as dividends above the threshold 
b (the name barrier is replaced by threshold, since the process goes on evolving above b). When 
7 = 1, the two-step premium policy reduces to a reflecting barrier policy. 

Threshold strategies are motivated by the fact that the optimal dividend distribution under 
the De Finetti objective with an extra constraint 

dD{t) < d dt, 

where 0 < d < c, is known to be consist in modifying the premium fromc to d above a constant 
barrier b - see for example [ATOQt IAM12| . Another motivation is provided in |SYY13j . 

For CramCT-Lundberg processes, ruin is sure or not depending on whether [BLPll] 


1 - 7 < d : = 


A^[Ci] 

c 


( 12 ) 
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or not, where A, c denote the rate of arrivals and premium, respectively. 

For threshold/refraction policies, some formulas expressed in terms of scale functions are 
provided in [KLlOt |Kypl3| . 

3. A third possibility is that the subsidiary pays ’’tax” - a proportion 7 of the premium income, 
whenever the subsidiary is at a running maximum- see for example |ARZ08l IAACI141IAI14] . 
(when 7 = 1 , this reduces again to paying dividends above a constant reflecting barrier). 

We examine now one possible efficiency criterion. 

3.1 Efficiency based on De Finetti constant dividend barrier and linear penalty 


As noted above, a critical step in defining efficiency is the choice of the value function. Our choice is 
based on the case of De Finetti constant barrier and linear penalty (with k proportional cost and K 
the fixed cost), for which the value function is known explicitly, see |APP09l (13)], [LRlOl lAPPlSj . 

Putting F{x) = Jq F{u)du,F{u) = k {Z{u) — k!{Q)W( u)) — KqW{u)^ , they find that the value 
function may be written as 


V{x) = F{x) + W{x)G{b), G{b) 


IF'( 6 ) 


1 - k {Z{b) - K'{0)W{b)) + KqW{b) 
W'{b) 


The optimality condition is obtained by differentiating the ’’barrier influence function” G{b) 
|APP15t Sec. 13.1]. Putting q := A(6) = c (w"^{b) — W{b)W"{b)J, we hnd 

H{b) := G'{b){W'{b)f = {k (k'(0)1F'(6) - qW{b))) + KqW'{h)) W'{h) 

-W”{h) (1 + k{K{Q)W{b) - 1 - qW(b)) + KqW(b)) 

= (Kq + k(l - p)) c (w'^{b) - 1F(6)1F"(6)) + kq {W"{b)W{b) - W\b)W{b)) + {k - 1)1F"(6) 
= k{{l- p)A{b) + q {W''{b)W{b) - W'{b)W{b)) + W"{b)) + KqA{b) - W”{b). (13) 

Now to compute 1F"(0), let K{s) = ^ ^ \ ._ A the expression used in [?, pg. 33] to 

^ l—XF{s)—S/s 

find IF'(O) = lims^oo Ar(s) = Following the same approach, we find 

1F"(0) = lim s{K{s) - W\0)) = i ((A + ^2 _ a/^(o)) . ( 14 ) 

5^CX> C \ / 


Using furthermore 1U(0) = 0, VU(0) = 1/c, cA(0) = A/c(0), we find: 


cG'{0){W'{0)f = k(^{l- p)Xfc{0) - qW'{<S) + clF"(0)) + Kf\fc{<S) - clF"(0) 
= k (a(A + g) - pXfciO)) + KfXfciO) - ((A + g)^ - A/c(0)) 


The optimality condition is: 


G'(0) <0^k< f{q,K) : = 


~,_, iX + ^^-XfciO)-KfXfciO) 

g + A(1 - mi/c(0)) 


(15) 


^F{u) = -KZ^°\u) + kZ^^\u), where the coefhcients of —K,k in F{x) are found by differentiating the second 
scale function | APP15 | Z{x,9) = (l — k{ 9) e ^^W{y)dy), 0 and 1 times respectively, with respect to 9, and 
then differentiating once more with respect to x. 
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Remark 5. Note that if K = 0, fc{0) = 0, this reduces to SLG readiness. 

Example 4. As a check, with exponential claims, the scale function is: 

where A± = and C,'^{q) = ^((?), C,~{q) are the largest and smallest roots of the 

polynomial c~^{'iIj{s) — q){s + p) = s"^ — s{X + q — p) — qp: 


Thus, 


C^{q) = 


q + X- p±\/{q + X-p)‘^+ ^qp 


1E"(0) = A+{C+{q)f - A_{Q-{q)f = MA+ -A.) + {X + q-p) (A+C+(g) - A_C{q)) 


and 


cVE"(0) = qp + {X + q - p){X + ^ = {X + - Xp, 

confirming m- 

Finally, with K = 0 and exponential claims of rate p, we recover JAPP 151 Lem 13.2 b]: 

^ ^ (A + - Xp ^ (A + qf - Xcp 

~ Xq Xq 

Lemma 1. With general claims admitting a first moment such that p < mifc{0) < 1, the function 
f{q, K),q>0 defined in ([I5D is increasing for any D < K < + rnJdo) ■ 

Proof: The function defined in the right hand side of (|15p is eventually increasing when q —?• oo. 
Its value in 0 is proportional to and hence negative under the hypothesis. The derivative 

with respect to q of of (fT^ is proportional to 

[z + l- /c(O)mi)^ - fci^fml + /c(0)mip"^ +Kfc{Qi) (/c(0)mi -l),z = q/X. 

Finally, assuming p < 1, we find three different conditions ensuring that the derivative is positive 
for any z > 0: 

K/mi < i_^^'7c(0) + FNho) /c(0)mi<l 
< K >0 fc{0)mi = 1 

K/mi > /c(0)mi > 1 

Moreover, in the case /c(0)ut,i > 1 if we have K/mi < then there are ri{K) < 

r 2 {K) such that the function q —)• f{q, K) defined in (IlSp is increasing on [0, ri], decreasing in [ri.r 2 ] 
and increasing to oo on [r 2 , oo). Because the image of g —>• f{q, K) is [0, oo) we can always find an 
appropriate killing rate, e.g. the smallest solution to k = f{q,K) (if k G r 2 ]), there are 

three such solutions, and we can choose the smallest one). 
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Theorem 1. With claims satisfying /c(0)mi = 1 (for example exponential), the readiness function 
is: 


f{Q) 


jX + q)^ - ^ ^ a 2mi - K ^ Ami - 1 

Xq A mi miq 


Three possibilities arise for subsidiaries. When K = 0, they are: 


1. The loading condition p < 1 is not satisfied, resulting in rejection. 

2. The loading condition p <1 is satisfied, and k < resulting in acceptance. 


3. 


The loading condition p < 1 is satisfied, and k > 

level £{ri), where r* = 5{k) — q, and S{k) is the solution of k 


, resulting in bailouts up to the 

_ (\+Sfi —Xcfi 

~ \S ■ 


Proof: The derivative of the right hand side of (fTOl) is proportional to <5^ + X{pc — X), f{q, K) is 
increasing whenever the loading condition p < 1 is satisfied, and therefore an exponential subsidiary 
not satisfying the efficiency condition (jl6p can always be rendered efficient by an appropriate killing, 
increasing its discount rate. 


4 Central branches as Sparre-Andersen processes 

Consider a central branch with SA subsidiaries that infuses capital into subsidiary i every time 
its surplus level drops below 0 for the j-th time. For simplicity, we assume from now on that the 
subsidiary is reset to zero. Given the challenges presented by multidimensional CB networks, we 
focus now on exponential or phase-type subsidiary claims, and the latter are general enough for 
our purposes. 

Remark 6. Taking Sfit) to be spectrally positive Levy perturbed compound processes SPLPCP 
\Fro08\.rLWS09[lZYY13^ (for example with a Brownian motion perturbation) poses often no problem. 
In this case, the subsidiaries may be kept nonnegative using minimal Skorohod regulation. Then: 

Xo{t) = uq + Cot - So, (17) 

/ 

5*0 = V' kilfit), Ii{t) = - inf{Xj(s), 0}, 

S<t 

i=l 

where the regulator process fi is the minimal process whose addition to Xi ensures that the sum is 
non-negative. 

With (/3, B) subsidiary claims, the bailout (time, size) pairs {Ak, Ck)k are IID random variables 
with joint distributions of the special form 

P[Ak G dt, Ck G dx) = d{t) 6'®'=^ bk dx dt, t, x G M+, 

where Bk := k~^B, and a{t) = {ai{t),..., ai{t),...) contains the densities of the ladder time joint 
with ruin in phase i. 

The CB is itself a Sparre-Andersen process with phase-type claims, exhibiting however a non¬ 
standard dependence (jl8j) : the arrival times of this process are the first-passage times for the 
subsidiaries and the jump sizes are given by the bail-out amounts (undershoots). 
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Remark 7. The classic Sparre-Andersen model with jumps of phase-type {I3,B) and independent 
inter-arrivals with density a{t) is a particular case of ([TH]) . obtained by taking densities of the form 

d{t) = a(t)P 

The essential difference in our case from the classic independent Sparre-Andersen process is that 
here the initial phase of a service (claim) period is decided at the bottom of the up-jump representing 
its inter-arrival time, and decides therefore also the size of the jump. 

Remark 8. Consider a CB network with several independent subsidiaries starting all at Ui = 
0, ,z = 1 ,...,/, having claims of phase-type = 1 ,...,/, and let p^^\t), R^^\t) denote the 

respective down ladder densities and survival functions. 

By conditioning, we find that the density of the the first bailout of the CB is 

P{Ak € dt, Ck S dx) = I R^^\t) I (t) e^^i * dxdt, t,x & M+ (18) 

*=i \j¥^i ) 

where a^'^\t) have Laplace transforms satisfying Kendall equations. 

In the case of one subsidiary, 1 = 1 the above equation gives the density for all bail-outs. When 
/ > 1 this is no longer the case since only one of the risk processes is controlled at the time of the 
first bailout. 

The Sparre-Andersen process representing the CB is considerably harder to analyze when I > 1. 
In the sequel, we analyze the CB process in a stylized case of hierarchical networks for I > 1 and, 
respectively, we provide approximations of many quantities of interest based on the CB process in 
the case 1 = 1. 

5 Linear networks of deterministic CB: reduction to one dimen¬ 
sion 

It turns out that as long as the CB in isolation is a deterministic drift with parameters uo,co, 
the ruin probability of the CB equals that of a subsidiary with modified parameters 

u[ = ui + uo/k, c'l = Cl + Co/A:, 

where k = ki, independent of the reset policy! 

The result is the same as if the CB transfers everything at the time 0+. This is also the case 
with several other problems involving a drift CB with no extra liabilities, which ends up liquidated 
totally. 

The pooled assets auxiliary process. Introduce the pooled assets combining Xq and the 
reflected processes Xi{t) = Xi{t) + lift) in such a way that the transfers and regulation cancel out: 

I 

X{t) = Xo{t) + Y,kiiXiit) + R{t)) 

i=l 

I 

= uq + CQt+ '^kiXi{f) 
i=\ 

I 

= u -\- ct — kiS^^^ (t) 
i=l 
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where we put 


1 I 

u = UQ + '^kiUi, c = cq + y^^kjCj. 
i=l i=l 

Remark 9. Note that with spectrally negative Levy subsidiaries X{t) is also a spectrally negative 
Levy process, while Xo(t) is a complicated superposition of SA processes. 

However, when 1 = 1, the ruin time of Xq and U coincide. Furthermore, the pooled reserves 
from the point of view of the subsidiary 

has the same law as the subsidiary with combined initial value and income rate! 

Remark 10. With several subsidiaries, 

T = inf{t : X{t) < 0}. (19) 

represents the ruin time if the subsidiaries may start helping each other at no cost, once the CB is 
ruined. 

These remarks yield the following: 

Theorem 2. Let Xq be a CB with deterministic drift and arbitrary structure subsidiaries. 

i) Assume 1 = 1 and put k = ki. Then, the ruin time tq of the MRN equals a.s. and in 
distribution the time r of the pooled process defined in dni), and equals furthermore the ruin time 
of the subsidiary with modified initial reserve UQ/k + ui and premium rate co/k + ci, Fi{x), namely 

^{u,c,t) = ^i{uQ/k + ui,CQ/k + ci,t), (20) 

independently of the reset policy! 

a) For I > 1, the time r is an upper bound for the ruin time tq of the CB: 

To < T. 

Hi) Statements similar to A) hold for any Cerber-Shiu objective, with or without dividends to 
the subsidiary, as long as the CB is a deterministic drift. 

Optimal allocation of total reserves = uq + ui and premium rate c+ = cq + ci. Note 
that: 

1. When /c > 1, uo = Co = 0, ui = X, Cl = c achieve the minimal ruin probability. 

2. For k = 1, the ruin probability is independent of the amount ui G [0,tt+], as well as of the 
amount of premium rate ci G [0, c+]. 

This optimization result hts the intuitively clear fact that with one subsidiairy and no expenses, 
it is optimal to take advantage of the hrst transfer without cost to transfer everything to the 
subsidiairy. 

The next corollary extends the previous result to a linear chain of CB, which can represent a 
chain of reinsurers, see [?]. 
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Corollary 1. Let Xq, Xj-i denote a linear chain of CB’s with deterministic drift. Assume 
Xi,i = 0,— 1 must pay proportional costs ki for bailing out Xj+i. 

Then, the probability of ruin of the MRN satisfies 


^{u,c,t) = 


Uq 


+ 


Ui 


+ ...ui, 


Co 


+ 


Cl 


'ko...ki-i ki...kj-i ko...kj_i ki...kj_i 

where ^j is the ruin probability of the last subsidiary in the chain. 


+ ■■■ci,t), 


Remark 11. One interesting feature of this result is that it does not require any assumption on 
the probabilistic structure of the subsidiary risk process. 

Another interesting feature is that similar reductions hold for other problems, as long as I = 1 
and the main branch is a deterministic drift (in the absence of subsidiaries). For example, one may 
add subsidiary dividends, ruin observed only at Poissonian times, Parisian ruin, etc. 

Example 5. Consider a CB with a spectrally two sided Levy subsidiary X{t), and ruin time 
To observed only at Poissonian times Tj,i = 1, ... of frequency rj. Recall M/i5L (2)J that the 
survival probability in this case is given by 


(pniu) = E4>{u + U), 

where U is the up factor of the Wiener-Hopf decomposition 


U = sup X{t) = U{S^) - inf X{t) 


Supposing furthermore that the subsidiary owes an exponential ruin severity penalty, one 
may use JAI151 (6)] that the survival probability in this case is given by 


E; 


U.T} 


g-gTo+0[/(To). < oo) =T 


e Eu+U 


g-gro+0t/(ro).^^ < OO 


E^-,T^+eD 


where D = info<s <£:,,is the down factor of the Wiener-Hopf decomposition, and ,T^ are 
the first times occurrences of U,D. 

Note that dividends over a harrier, and/or Parisian ruin may be added as well, using \ All 51 
Prop. 4], o,nd \AI1,% Prop. 4]- 

If furthermore the Levy subsidiary is spectrally negative, the joint Laplace transform of the 
Poissonian observed ruin time and of the ruin severity is more explicit jAIZlfl (14)] 

< oo] = 

{z(u, 0) - z(u, (/ - 

Finally, the ruin probability of the CB will be given by the same formulas, with modified initial 
capital u and premium c. 


Further questions of interest, not considered here, are: 

Problem 1. Determining the optimal allocation of the initial total capital 
total income rate c+ = J/j^gCi, when I >2, for example under red time objectives. 

Problem 2. Another important question is how to dynamically manage the subsidiaries, i.e. 

(j\ 

how to set the reset levels x) optimally, possibly adaptively, depending on the evolution of the 
information available. 


Problem 3. Study CB branches with an extra Levy spectrally negative perturbation: for example 
with exponentially distributed claims, or with Brownian motion (BM). The problem is particularly 
important, since the BM could model an approximation of other subsidiaries, and is considered in 
the next section. 
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6 SNMAP approximations of non-deterministic CB 

In the reminder of the paper we focus on the case of a non-deterministic CB with one subsidiary. 
We have seen in the previous sections that the risk process of the CB is itself a Sparre-Andersen 
process in which the inter-arrival times are the hrst passage times and respectively the bailout 
amounts for the subsidiary. We provide here a methodology to arrive at an SNMAP approximation 
the CB risk process, with the goal to provide an input to the available scale methodology toolbox 
l |KP08l [Ivain IIP121 UvalSj i which can be used to simultaneously solve a wide variety of insurance 
problems for this CB risk process. We exemplify in this section with the computation of the optimal 
dividend barrier for the CB itself, and implicitly for the whole network. 


6.1 Two point Fade approximations for the downward ladder time of the Cramer- 
Lundberg process with exponential claims 

In this subsection we consider the particular case when the subsidiary risk process follows a Cramer- 
Lundberg process with exponential claims and we are interested in its downward ladder time. 

Our approach is based on the idea that approximating the excursions of a process ensures 
approximating the process, and in particular various functionals of the process |LSZ13l lYanl3] . 

With phase-type jumps, one would need to provide bivariate matrix-exponential approx¬ 
imations for the joint law of the downward ladder time and height of a SA process. 

In the case of exponential claims, the density of the downward ladder time may be expressed 
as a hyper geometric function: 

p{t) = p{c^it)cn, p{t) = pe~^^+P'>\Fi{2,pt^). (21) 


However, what we need is a phase-type approximation of this. This topic has already been 
considered in |AW88j (at order two), as one of many possible methods for approximating the M/M/1 
busy period density. 

We recall now some basic facts on this case |AM15| . 

Lemma 2. a) With Poisson arrivals of rate A and exponential elaims of rate p, the Laplace 
transform of the downward ladder time density satisfies a quadratic equation 


Pq = a{q + cp- cpps 


_ A 

X + q + cp - cpps 


P 

1 I X I j P 

r P 0 -\- p — ps 



Q 

cp 


with solution 



(^6 + p+1- -^^{5 + p+iy - . 


( 22 ) 


b) The Laplace transform (12^ may be computed iteratively by the continued fraction expansion 


Ps = 


1 + 6 + p- 


l-h5+p—s... 


The convergents pWn = 1,2,3... with s constant satisfy p\ ^ \Neu6(k (75)] where 

Pn{x) are Chebyshev polynomials \Neu6(k (77)J. When s = 0, the first three are 


p{6 + p + l) 


p + 2p6 -\- 25 -\- p^ -\- p ^ 


5 + p + F 5“^ + 2p5 + 25 + p^ + p+F (<5 -I- p -I-1) (d^ -I- 2p5 + 25 + p"^ + 1)' 
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Decomposing in partial fraction the third convergent yields an order three rational approxima¬ 
tion of the Laplace transform of the ladder time: 

p {5 + p+lf-p 


PS 

P 

2 


5 + p + I {5 + p + ly - 2p 
1 1/2 






1/2 




(23) 


5 -\- p \ 5 -\- p \ \/2p 5 -\- p \ — y/2p ^ 

Inverting the Laplace transform yields a hyper exponential density approximation: 

m « + e+v^*)) := 

i=0 

_ i \/2p)t 

2 ^4 ^4 ’ 

where ao = 2 {i+p) = 4 (i+p^,y 2 p) nonnegative for any p S [0,oo). Furthermore, 

any p in this range. 

may he obtained faetoring 

6)-^: 

6 + p+1 


iff P < Ij providing us thus with a valid approximations for any p in this range. 

Remark 12. A further simplification of the Laplace transform 
5 + p + 1 and ehanging variables a = p{l + p + 6)~‘^: 


PS = Ps = 


P 


1-.1-4 


P 


{5 + p + 1)2 


1 - VI -4a 


(24) 

h -|- p -|- 1 2a 

The second factor put thus in evidence is the generating funetion of the famous Catalan numbers 
1 - VI-4a °° 


2a 


(2k\ 

= 1 + a + 2a2 + ha^ + 14a^ + ... 


k=0 


k + 1 


and a continued fraction (cf) representation 

1 - yi - 4a 
2a 


1 - 


(25) 


1 - 


1 + ... 


may be found for example in \CPV^ 0^ (7.7.5)]. H The lowest order approximations dll]) are 

1 1 — a 1 — 2a 1 — 3a-|-a2 1 — 4a-|- 3a^ 

1 — a’l — 2a’l — 3a-|-a2’l — 4a-|- 3a2 ’ 1 — 5a -|- Qa? — af" ’ 

Lemma 3. The rational convergents Rn of the continued fraction (l25]) increase towards Va € 

(0,1/4). 

Proof: This is immediate by the positivity of a. 


Remark 13. Alternatively, we may use two point Fade approximations which ensure also the 
equality of the derivatives around 0. \AW8^ Sec 3] provide an in-depth numerical comparison of 
several hyper-exponential approximations of order two, and find that fitting the derivatives yields 
excellent results around 0, while fitting the moments is less satisfactory, since better results may be 
obtained with asymptotic approximations. 

^ As well known, the Fade approximations obtained by truncating continued fractions have good properties, like 
larger domains of convergence than the corresponding power series. 
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Let us invert now the lowest order two-point Fade approximation of p{6) which ensures also the 
condition = p, p £ (0,1]: 


p ((5^ + 2/9(5 -\- ‘15 -\- p 


P 


((5 + p + 1) (,52 + 2/9(5 + 2(5 + 1) p + 2 + P + 1 <^2 + 2(5 (p + 1) + 1 


P + 1 




S + p+1 


(26) 


P 


P + 1 




1 


P + 2 y(5-|-p-|-l Ai — A 


- (p+ 1) p-M - A; 


(5 + Ai 


(5 -|- A 2 


^ , Ap2 = P + 1 + Vp(p + 2) 


This yields a density approximation: 


pit) ~ 




(27) 


i=0 


where , ai = 


p 


i-(p+i)/Ai _ p^^^(Vp(p+2)+i) 


p+2)'-i p+2 A 1 -A 2 2 ( p + 2 ) 3/2 (^p+l+^p(p+2)y ’ 


0:2 = 


_ p (p+ 1)/A2-1 _ P^^^(\/p(P+2)-l) 


p+2 Ai—A 2 


and ai = p. 


2(p+2)3/2 ^p+1—p(p+2)^ 


Note that ai are nonnegative for any p e [0, 00 ), 



Figure 1: Two order 3 matrix exponential approximations (I23p . (j27|) of pit) = ipit] 0) with A = p = 
l,p = I = c~^. The Bessel law (blue) is in between the two point Fade approximation (in red), 
which is tight both at 0 and 00 , and the classic continued fraction Fade approximation (in yellow), 
which is better near 0 


6.2 Fade based SNMAP approximations for A(t) when subsidiary claims are 
exponential 

Having established in the previous subsection convenient approximations for the subsidiary risk 
process, we now approximate the CB risk process. 

Note that even though finite time ruin probabilities have an explicit Bessel density with expo¬ 
nential claims, we will replace them by matrix exponential approximations, since this allows solving 
network problems by the SNMAF methodology. 

After applying the order three approximation (|27p to the subsidiary’s ladder time, the central 
branch becomes a MAF with three states, with transition rates Qij = Xiaj accompanied by expo¬ 
nential jumps of rate p/k translated by K (we could include here phase-type jumps to the CB, 
and fixed costs, since these pose no problem to the MAF methodology). 
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Figure 2: Two order 3 matrix exponential approximations (l23|) . ( 1771 ) of 0) with X = = 1, p = 

^ = c~^. The fit is better as p gets further from 1, and could be improved by using a combination 
of the two approximations, of order 5. 

When K = Q, f(s) = and we find from the general formula that the symbol of the 

approximated CB is 

]K(s) = diag{cos — A*) + Aa/(s) (28) 

^ cos - Xo + Xoao^^ ^ 

= cos - Xi + Xiai^^ 

y A2ao40fc cos-X2 + X2a2j^ j 

To apply the scale based MAP methodology, it is convenient to transform this MAP with 
exponential jumps into a continuous MMBM: 


/ cqs — Ao 

0 

0 

Aotto 

Aotti 

Aoa2 \ 

0 

cqs — Ai 

0 

Aiao 

Aicri 

Aia2 

0 

0 

cqs — X2 

A2Q;o 

A2Q:i 

A2a2 

p/k 

0 

0 

—s — p/k 

0 

0 

0 

pjk 

0 

0 

— s — p/k 

0 

V 0 

0 

p/k 

0 

0 

— s — p/k / 


(29) 


After obtaining an SNMAP approximation for X{t), we may solve approximatively various 
problems related to this process, using the package |Iva m- 

We illustrate this by finding in Figure [3] the optimal dividend barrier. It is beyond the scope of 
this article to investigate further this problem in particular (which deserves a separate treatment), 
as our main goal was not the solution of the the problem per se but providing an adequate input 
to the toolbox that solves this type of problems. 


7 Spectrally-negative Markov additive processes (SNMAP)- a new 
framework for risk modelling 

Having demonstrated the power of the approximation based methodology, we dedicate the remain¬ 
der of the paper to set bases of future research directions that build on the extraordinary potential 
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Figure 3: Dividends as function of barrier, starting from ^0 = ^ 1 = 0, with A = // = A: = 1, 
Cl = 4, Co = 24, obtained by the scale matrix methodology [Iva m. The maximizing barrier is ~ .17 

in insurance of the class of spectrally-negative MAPs (SNMAP), that is, MAPs which have negative 
jumps only. 

Informally, a MAP is a bivariate process where J{t) is a Markov chain (MC) rep¬ 

resenting an exogenous background process, and X{t) is a so-called additive component modulated 
by J{t) (nonetheless we often say MAP meaning X{t)). 

A MAP is a generalization of a Levy process in the sense that X{t) has stationary and inde¬ 
pendent increments conditioned on the state of the modulating process J{t). 

Definition 2. A bivariate process {X{t), J{t)) is called MAP if, given {Jt = i}, the shifted process 
{Xx+t — Xt] Jx+t) is independent from {X(t), J(t)), 0 <t<T and has the same law as {X{t) — 
A(0); J{t)) given {J(0) = i} for all i and T > 0. 

A spectrally-negative MAP is a MAP whose additive component can have only negative jumps. 

Remark 14. Furthermore, one can replace the deterministic T in the above definition by a stopping 
time. The resulting property is called the strong Markov property for MAPs. 

It is common to assume that ,J{t) is an irreducible MC with a finite state space ff = {1,... ,n}, 
which we do throughout this work. It can be shown that X(t) evolves as some Levy process XI 
while J{t) = i, and in addition a transition of J{t) from i to j may trigger a jump of X{t) distributed 
as Uij, where J{t) and all the components in the construction are assumed to be independent. This 
construction presents an alternative often-used definition of a MAP. See |Asm03[ Ch. XI] for further 
information on MAPs. 


Remark 15. In the presence of a ’’modulating environment” J, one key idea is to consider matrices 
of expectations and probabilities, conditioned on some phase i at start and joint with some phase j 
at crossing the level of interest. These will be denoted E = E[...; J(t)],P = P[...; J(t)]. 


The Laplace exponent. The law of a spectrally-negative MAP is characterized by a certain 
matrix-valued function IK(0), defined by 


E[e®^W; J(t)] = 
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where the (i,j)-th element of the matrix on the left is given by J{t) = j\J{0) = i) H . 

The n X n matrix ]K(0) is given by 


' \qijU,,ie), ifi/i, 

where qij are the elements of the transition rate matrix Q of J(t), Ki{0) = logEe^^®^^^ is the Laplace 
exponent of the Levy process XI, and Uij{6) = In matrix form 

K(0) = Diag{Ki{9)) + Q0U{9) 


where 0 denotes Hadamard product. 

Remark 16. The Levy processes may have killing parameters qi = —Ki(0) > 0. When Ki(0) = 0, 
it follows that ]1C(0) = Q is the transition rate matrix of J. 

Remark 17. An alternative notation in use when a fixed nonzero killing q is involved is to exclude 
it from the symbol matrix IC.(0), and write Kq{9) = K{9) — ql when killing is present. We will 
follow, for the sake of readability, the convention of including a fixed killing argument q in the 
symbol, scale, etc and just write K, VL,... instead of Kg, Wg,.... 

If an additional secondary killing g appears in certain circumstances, its presence must be indi¬ 
cated. The notations ]Kj.,ILr,... corresponds then to VFr+g,... in the complete notation. 

Example 6. A Markov-modulated Cramer-Lundberg process with tax is retrieved from a spectrally- 
negative MAP {X{t), J{t)) by putting Uij = 0 and XI = Cjt — where Ni{t) is a Poisson 

process of intensity A, and Ck are iid positive random variables. Hence Gij{9) = 1 and Ki{9) = 
Ci 0 - Ai(l -Ee-^^i). 

Remark 18. With phase-type jumps, it may be more convenient to work with the symbol of the 
’’embedded process”, in which jumps have been replaced by downward drifts M.sm.9,5l \ AA1(^ . 

In applications, we are interested in ’’regulated” versions of X(t): 

X[o(t) = X{t) + Lo{t), (t) = X{t) - Rb{t), 

X^t) = = Xit) + L{t) - ^R{t), 7 < 1, 

where Lo(t) = —(X(t) AO), 2£.{t) = info<s<tis the minimal regulator constraining X{t) to be 
nonnegative, 

Rb{t) = {X{t) - b), , X{t) := sup X{t) 

0<s<t 

is the minimal regulator constraining X{t) to be smaller than b. The process is said to 

be reflected at 0 and refracted with coefficient 7 [AI14) (1)]. The minimal regulators L{t),R{t) 
have points of increase contained in {t > 0 ; T(t) = 0} and {t > 0 : Y{t) = Y{t) V b} respectively, 
and their explicit formulas are more complicated in this case; for a recursive construction, see for 
example [APPOTl pg 165]. The process will be denoted by X^{t). 

Remark 19. The process 7 l?(t) intervenes here as cumulative dividends paid to some beneficiary, 
and using 7 < 1 allows X'^{t) to continue above the upper ’’threshold” h (with a modified drift). 

^This is the analogue of the ’’Laplace exponent/symbol/cumulant generating function” of a spectrally-negative 
Levy process. 


21 











Smooth exit upwards and the associated generator. One key idea in the study of 
SNMAPs is considering the phase process ’’while progressing upwards” J+ , which is itself a MC as 

_ ‘x 

a function of x > 0, by Definition [2] with T = r+. The transition rate matrix G is a (sub)generator, 
the matrix 


P[J,+] = (30) 

contains the probabilities of hitting x in various phases, and 

contains the probabilities of ever hitting x, starting from all initial states. 

The matrix G is a right solution of the equation IC(—G) = 0 |Bre08] . 

Remark 20. In the Levy literature, G is a scalar denoted by —which solves the Cramer-Lundberg 
equation k(<P) = 0, > 0. 

If the MC J^+ is transient, then its life time coincides with the overall supremum of X{t), and 
has thus a PH distribution characterized by the matrix G. 

The two-sided smooth exit problem and the scale function. Cf. HEE], the 

solution of the smooth two-sided SNMAP exit problem has a multiplicative form (well-known in 
the case of Levy processes): 

Theorem 3. There exists a unique continuous function W : [0,oo) — )■ such that 

e-^^W{x)dx = K(0)-^ 

for all sufficiently large 6 > 0, W(x) is invertible for all x > 0, and 

<TQ]J{T^)] = W{x)W{b)~^forallb>0,xe[0,b]. (31) 

The decomposition 



W{x) = e-^^H{x), 


where H (x) is the matrix of expected occupation times at 0 up to the first passage over x, provides 
a probabilistic interpretation of the scale matrix. 


Remark 21. Informally, W{x) is the solution of 
generator of the process X{t). 


Giy(x) = 0, 
W{x) = 0, 


X > 0 

X < o’ 


where G is the Markovian 


Remark 22. A matrix expansion for the two-sided exit by jump probability with phase-type claims 
and inter-arrivals was already obtained in \APU0A, (13)-(14)] via an elementary ODE approach, 
later extended in IJP081 (26)-(28)]. However, the intuition and rigorous proof that everything 
could be expressed in terms of one scale matrix were furnished only later, in the seminal papers 
\KP08\. \IP12I] . respectively. 
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Example 7. A Sparre Andersen renewal risk process with inter-arrival times of phase-type {a, A), 
is an example of a spectrally-negative MAP {X{t),J{t)). Here the MC J{t) lives on n states, and 
J(0) is distributed according to a. For t > 0, J{t) makes a jump from i to j without causing a jump 
of X{t) with rate Oij; it makes a jump from i to j and triggers a jump —Ck with rate aiaj. Hence 
for i j it holds that Fij{9) = aij + aiaj¥,e~^^^, that is qij = Oij + aiUj and Uij is an appropriate 
mixture of 0 and —Ci. Then Fii{9) = an + aiaiKe~^^^ + c9, because qn = — dij = 
and Ki{9) = c9 — aiai{l — which corresponds to a compound Poisson process with intensity 

aiai, jumps distributed as —Ci and drift c. In matrix notation we have 

K(0) = ^ + aaEe-^^i + c9ln, 


where is an n x n identity matrix. 

The second scale function [APP151, IIP12| is defined by 

Z(x,9) = e^^ (^I - e-^yW{y)dyK{9)^ . 

Alternatively, defining the Dickson-Hipp transform 

roo / rx 

W^(9) = J e-^m{x + y)dy = i^{9)-^ - J e-^yW{y)dy 

where the equality holds for large enough, it holds that 

Z{x, 9) = 1T,,(0)IC(0) 4^ W,,(9) = Z(x, 9)K(9)-\ 


(32) 


i.e the second scale function Z(x, 9) coincides up to a constant matrix with the Laplace transform 
of the shifted scale function (the ’’normalization” ensures that Z(O,0) = I). 

Another way to characterize the second scale function is via its Laplace transform: 

Z{s,9) = {9- s)-^(IC(s)-^ - K(6')-^)IC(0) = {9 - s)-^(K(s)-^IC(0) - I). 


In the Levy case, this becomes 


Z{s,9) = 


k{s) — k{9) 


{s — 9)k,{s) 

The second scale function with two-step killing is |AI13j . |AIZ14l (12)]: 

< To"; j(r+)] = Zix,^r)Zib,^rr^yb > 0, X G [0,b], 


(33) 


in the sense that it replaces the first scale function in the smooth two-sided exit problem with extra 
killing below 0. 

Remark 23. The function Z{x,9) appeared first in 2008 in \APP15j . in the context of finding 
smooth Gerher-Shiu functions associated to an exponential payoff e^^. Subsequently, its key role 
in many other first passage problems was revealed in UP10j and fAIZlj^ . 

Killing when the draw-down exceeds a given value a. Consider a process X^\t) = 
X{t) — Rb{t) = X{t) — {X{t) — 6)+ regulated and started at b, resulting in 

a'’] (t) = X{t) + b - X{t) = 6 - y°(t). 
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where 


Y^{t) = x + X{t)-X{t), Y^{0)=x 

is called ” drawdown/reflection from the running maximum” (starting from x) |Kypl4 


pg. 248], 


The regulator Rb{t) = {X{t) — 6)+ can also be seen as the total amount of dividends paid until 
ruin in a Levy model with the barrier dividend strategy, where the initial capital and the barrier 
are both placed at the level b. 

Let us kill now X{t) (send it to some absorbing state) at the stopping time 


ta = inf{t > 0 : Tt > a}, 


i.e. at the first time when the height of an excursion from the maximum exceeds a > 0, or, 
equivalently, at the first time when the regulated process X^^ (t) started from b > 0 drops by more 
than b — a. 

Using the strong Markov property for MAPs, we see that also in the presence of killing at ta, the 
environment phase observed while evolving upwards J + ,x > 0 is still a MC, with some transition 

'^x 

rate matrix L(a), so that P[r 3 / < ta, J^+] = 

It was shown in |IP12[ 4], generalizing a well-known excursion theory relation for Levy processes, 
that for a > 0 the right and left derivatives IU'(a+) and W'_{a) exist and 

L(a) = -W'{a+)W{a)-\ L(a-) = -W'_{a)W{a)-\ P,[r+ < 4; J^+] = (34) 

Remark 24. For SNMAP’s, the life time of the transient MC characterized by L{a) is a matrix 
exponential law ill). In the case of a Levy process, this reduces to an exponential random variable 
with rate W[oj^)/W{a) dill) . 


7.1 Resolvents 

The q-resolvent measure, for any Borel set B G [a, b], may be expressed in terms of the scale function 

-Wg{x - a)Wg{b - y) 


E, 


Ta At/ 


-qt-^ 


{Xt&BY 


dt = 


1 


{y&B} 


Wq{b-a) 


- Wg{x - y) 


dy a <y <b-, 


see Theorem 8.7 of |Kypl4|. Letting a —)• — oo yields 


Eo. 


e '^^l{Xt&A}dt 1 = 


) = ^ -y)- Wg{x - y)] dy. 


(35) 


(36) 


Resolvent of doubly reflected Levy processes As shown in l |Pis03l Thm. 1]), 


a version 

of the g-potential measure of a doubly reflected Levy process with upper barrier b is Ul^{x,dy) = 
/o°° ^~'^^^x{X^{t) G dy) of X^ is given by 

Ug{x,dy) = vP{x,b)5b{dy) + vf{x,y)dy 

where 5b is the point-mass at b and 

iB{x,y) = Zq{x)Zg{b)-^Wa{b - y) - Wq{x - y), x,y € [0,b],y / b 

y = b. 


” q\ 

d'^ix,!) = Zq{x)Z'{b)-^Wq{0), 


(37) 


'qi-^/^q\ 

The next section illustrates further the fact that the answers to a large variety of first passage 
problems for SNMAP’s may ergonomically be expressed in terms of the matrices K, G, H, W and 

Z. 


24 














8 A compendium of first passage formulas for SNMAP and SNLevy 
processes 

In the research avenue we envision, diverse problems in insurance and in particular multidimensional 
problems which so far have defied exact analysis can be solved using approximations by SNMAP 
and ’’the eight pillars of one-dimensional hrst passage problems for SNMAP”. We assemble here 
these key formulas, which have appeared in various prior works |APP07l IAIZ141IAI141 lAIlSj for 
the Levy case and |IP12l llvall| for the SNMAP case. 


8.1 Homogeneous problems 

1. The De Finetti expected discounted dividends over a barrier 6, satisfy a relation 
similar to m-- 


Ex 



e-'i^dRb{s) 


W,(x)W'(6)-\ 


(38) 


when Tq is the classic ruin time. 

With Poissonian observed ruin, the expected discounted dividends over a barrier b become 


IATZ141 (27)] 


Ex 



e-^^dRb{s) 


Zg(x,4>g+^)Z'(6, 4>g+^) G [0,6]. 


(39) 


In the case of continuous observations r ^ oo, this expression reduces to the previous result 

dSI). 

The law of the total dividends until ruin R[,[tq) is known as well; it is a matrix exponential 
distribution generated by L(6) = —W'(b)W(b)~^ - see (09]). 


2. Expected discounted dividends until the total bail-outs of a reflected process 
surpass an exponential variable Sq satisfy (see |AI14( (15)] for the particular case b = x) 


F = IE, 


/ e ‘^^MLo{s)<eg]dRb{s) 
J [0,cxd] 


= z,{x,e)z'{b,e) 


-1 


(40) 


When 0 = 0, this yields |APP07[ (4.3)], and when 0 —)• oo, this recovers (f38l) . 

3. Bail-outs of a reflected process, until the first dividend. Let -^[o(t) denote a SNMAP 
process reflected at 0, let L^it) = —(0 A A(t)) denote its regulator at 0, so that 2f[o(t) = 
X{t) + and let E^ denote expectation for the process reflected at 0. Then, the joint 

law of the time until a process reflected at the infimum climbs to an upper level 6 and of the 
bail-out (regulation) is: 


E°[e-^T+-e^o(r+)] 


Zg{x,9)Zg{b,6) ^ 0 < oo 

E[r6+ < Tq-] = Wgix)Wg{b)-^ 0 = OO ’ 


where Zg{x,6) is the second scale function |IP12[ Thm 2]. 


( 41 ) 
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Furthermore, [AI14[ (3)] show in the Levy case that a power relation holds when replacing 
X[o(t) by a refracted process started at b. The proof uses the probabilistic interpre¬ 

tation 


where Kg is the hrst moment when the total bail-out exceeds an independent exponential rv. 
£e- 

Finally, |AIZ14[ (22)] extend to the case when is replaced by 

8.2 Non-homogeneous problems 

1. The severity of ruin before seeing a barrier b (two-sided exit). Applying the previous 
result gU) and the two-sided exit formula ()3ip . one finds that [IP 121 Cor 3]: 

J(r+)] =Z(x,0)Z(6,0)-i = 

< To"; J(r+)] +E,[e^^("o-);r- < r+] E[][e-®^oK+); J(r+)] 

= VF(x)lF(6)-^ + E,[e''^("o“);r- < r+]Z(6,0)-\ 

We may solve now for the joint Laplace transform of the hrst passage time of 0, and the 
undershoot 

E, (e-''"o-+e^(A ); r- < r+) = Z,(x, 6) - Wq{x)Wq{b)-^Zq{b, 9) (42) 

In the SNMAP case, denote by 

H = lim H{x) 

X^OO 

the matrix of total expected occupation times at 0. Assume that either Q1 ^ 0, or k'(0) ^ 0. 
Letting 6 —)• oo above, and using [IP 121 Cor 4] 

lim W{b)-^Z{b, 9) = {R + ei)-^K{e) := K(0), (43) 

b—^oo 

where R = H~^GH is a left solution of the equation K(—i?) = 0, we hnd the severity of 
ruin (non-smooth one-sided exit): 

Va:,0>O, E,,[e®^("'o");ro" < oo] = Z{x,e)-W{x){R +ei)-^K{e) 

= Z{x, 9) - W{x)H-^{G + 9I)-^H^{9). 

Remark 25. Note ^hat if 9 is a zero of det{K{9)), then it is a zero of det{91 + G) and 
det{91 + R). Thus, K{9) should he interpreted in a limiting sense. 

In the Levy case, (1431) reduces to limb^ooW{b)~^Z{b,9) = yielding the classic [AIZ141 

( 7 )] 

E.[e-^"o-+^^Po-)] = z,ix,9) - W,ix)'^^^. (44) 
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Remark 26. When we apply the above with 6 = this becomes: 


E.r e 


-ITq +‘^rX _ 


( 1 + (g - r) / e-‘^^^Wg{z)dz - 


‘0 

rx 


0 ; Tg < OO = Zq{x, ^>r) - Wq{x) 


r — q 






^q - 


'•Wq{x) 


( 45 ) 


2. The severity of ruin at a Parisian ruin time, before seeing a barrier b. When Tq 
is replaced by a Parisian ruin time Tr, the severity of ruin satisfies in the Levy case IAIZ14[ 
(15)], [BPP^ (1.13)1 


E., 


^-qr,+eX{Tr)^^^ < r+ 


= I Z,{x, 0) - 0) 1 + ’■ - «(<>))'*■ (46) 


I 

When r ^ oo, this recovers (I42p . and when 6 —)■ oo it yields [AIZ141 (14)] 

Ex[e®'^^^o );ro" < oo] = (^Z{x, 9) - Z{x, 4»,.) ) r{q + r - ^{6))-^, 


by using 


Zq{b,^r+qr^Zq{b,9) 


^ {^r+q - ^q)K{0) 

b^oo r{9 — dig) 


(47) 


3. The joint law of dividends over an upper barrier, of the ruin time, and of the 
severity of ruin. Let denote the law of X reflected from above at b, and let Rb{t) denote 
the regulator. Prom |IP12l Thm 6], we find 

E6[g-»?ii6(ro-)+0X(r-);^- <^^] 

= Z{x, 9) - Wix) {W'{b+) + dW{b)) (Z'(x, 9) + dZ{b, 9)) (48) 

Remark 27. As a check, decompose in the two cases Tq < , Tq > rj^: 

E6[g-4?i?6Po-)+ex(r-).^- ^ ^ _ w{x)W{b)-^Z{b,9) + 

W{x)W{b)-^ (^Z{b, 9) - W{b) {W'^{h) + dW{h))~^ {Z'{x, 9) + dZ{h, 0))) 

= Z{x,9) - W{x) {W'^{b) + dW{h))~^ {Z'{x,9) + dZ{h,9)) 

With 0 = 0, X = 6, this yields the Laplace transform of the total dividends before 
ruin 

^\[e-^RArq).TQ < £q] = Z{b) + {W’{b+)W{h)-^ +my^ {-Z’{h)-dZ{h)) = 
{W'{h+)W{h)-^ + diy^ (W'{h+)W{hy^Z{h) - Z'{b)) = 

{W'{h+)W{hy^ + [{W'{h+)W{by^ - {W'{h+)W{by^W{b) - W{b)) K(0)] , 

^Note that in |LRZ131 ILRZ14I IBPPR |. the function Zq{x,^q+r) is denoted by Hq^r{x) = 

g'Sq+ra: ^ Wg(y)dy') 
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which is the Laplace transform of a matrix exponential law. For Markov modulated Levy 
processes without discounting, 5i = 0,Vi ]K(0) = 0, and this reduces further to 


<oo] = {W'{h+)W{h)-^+M) ^ W{b+)W{b)-^. (49) 

Thus, Ri){tq) has a matrix exponential distribution generated by L(6) = —W'j^{b)W {b)~^. 

When ?? = 0 and with equal discounting 5i = <5, Vi, (I48p yields the severity of ruin for a 
regulated process: 

= Z{x,e) - W{x)W'{b+)-^Z'{b,e). (50) 

This formula is also called the dividends-penalty identity [GLY06] . 

If furthermore 0 = 0, this yields the Laplace transform of the ruin time for a De 
Finetti regulated process, due in the Levy case to |AKP04] : 

Px[to" < £,] = Z{x) - W{x)W'{b+)-^Z'{h). (51) 

A further generalization in the Levy case with Poissonian observed ruin time Tg” is [AIZ141 
(23)]: 


)+eX{T^ ).J,- ^ ^ 

(Z(x, 9) - Z{x, $^)(($^ + V)Z(6, - rWib))-\Z'{x, 9) + dZ{b, 9))) . 


4. The Gerber-Shiu function common to the two-sided absorbed and reflected sever¬ 
ity of ruin problems. Note the similarity between the equations ()42p . ()50p for the Laplace 
transforms of the severity of ruin before seeing b, and with reflection at b, which involve the 
same functions VF(x), Z{x, 9). 

|APP15] show in the SNLevy case that this continues to be the case for any pay-off w which 
is ’’admissible” (satisfies certain integrability condition). 

Proposition 1. Given a < b < oo, x € (a, 5), and an admissible pay-off w : (—oo,a] ^ M, 
there exists a unique ’’smooth GS function” so that the following hold: 


V^{x) = E, (Y|,,b|(r-)) = F^{x -a)- W^’^^x - a) 


VM = Ex (A|„,fe](T-)))l = F^{x -a)- W^^\x - a) 


Fwjb - a) 
lT(9)(6-a)’ 
Kib-a) 


lF('')'(6-a) 


(52) 


Stated informally, this amounts to the fact that both these problems admit decompositions 
involving an identical ’’non homogeneous solution” 

The smooth GS function 


°° gk 

Z{x,9) = '^^Zk{x) 

k=0 

corresponding to w{x) = e^’’’’,x < 0 may be used as a generating function for finding GS 
functions Zk{x) associated to power payoffs w{x) = x^. Especially interesting are the cases 
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A: = 0,1 which intervene in problems with linear bailont costs. Using Z'{x,9) = 6Z — 
lU(x)]K(0), we find: 

Zo(x) := Z{x, 0) = Z{x), Zi(x) := = Z{x) - lU(x)K'(0). (53) 

The derivatives with respect to x are : 

Z'(x) = qW^{x), Z[{x) = = Z{x) - lU(x)K'(0). 


5. Further results when first passage is only monitored at Poisson times. Elegant 
generalizations of the classic first passage resnlts when observations are made at Poisson 
times with rate r are given in |AI13[ IAI141IAIZ141IAI15] (in the Levy case mostly): 

(a) Severity of ruin with additional killing above the barrier. Cf. [AIZ141 (18), 
Rem.2], in the Levy case it holds that: 


); To" < T+] = Tq" < oo] (54) 


When r ^ oo, this recovers (1421) . and when 6 —)• oo it recovers ()44p by (I47p . 

(b) The solntion of the non-smooth two-sided exit problem, when both exit times are mon¬ 
itored at Poisson times with eqnal frequency, is [AIZ141 (19)]: 


< T+j = Z{x,9) - Z(x,4>,) 


Wb[^r,9] 

Wb[^r,^r] 


K,{9){rl — k{9)) 


-1 


where f[x,y] denotes a Newton divided difference. 
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